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We develop a means of simulating the evolution and measurement of a multipartite quantum state 
under discrete or continuous evolution using another quantum system with states and operators 
lying in a real Hilbert space. This extends previous results which were unable to simulate local 
evolution and measurements with local operators and was limited to discrete evolution. We also 
detail applications to Bell inequalities and self-testing of quantum apparatus. 
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I. INTRODUCTION 

Coherence and interferometry are central to quantum 
physics, hence the fundamental importance of concepts 
like global and relative phases, single and multi-particle 
interferences, etc.. Phases are usually described by a 
complex number with norm one and the use of com- 
plex numbers seems intimately connected to the heart 
of quantum physics: the complex Hilbert space and the 
imaginary unit i in the Schrodinger equation. Still, we 
all know that the complex field C is isomorphic to a two 
dimensional real plane, and hence it is possible, at least 
formally, to compute all quantum predictions using only 
real numbers. It is even well known that it is possible 
to simulate unitary evolution by a quantum system re- 
stricted to real amplitudes and matrix entries by doubling 
the dimension of the Hilbert spaceQ. However, this sim- 
ulation breaks down for multipartite systems. 

Consider the case of two particles, first in a product 
state 4>®'^^ and let's add a phase. This phase can equiv- 
alently be attributed to either of the two particles: 

So far so good, a is only an (irrelevant) global phase, but 
this is no longer the case if the particles are in an entan- 
gled state. In such a case it is no longer clear to which 
of the two complex Hilbert spaces, the one describing 
the first particle or the second one, should be doubled in 
order to replace the imaginary part by additional real di- 
mensions. It is even not obvious whether such a doubling 
could "work locally" . By this we mean the following: as- 
sume the particles, two or possibly more, are spatially 
separated from each other. In such a case one may ask 
whether one can double some or all the Hilbert spaces 
and let each party (each party holds one and only one 
of the particles) manipulate the phases independently of 



each other, including in the case of arbitrary entangle- 
ment. 

Another situation in which the real simulation seems 
to break down is in the case of continuous time evolu- 
tion. It is generally understood that over time a state 
will evolve according to the Schrodinger equation, pick- 
ing up a complex phase due to the i in the exponent. It 
would seem that, although discrete time evolution may 
work with only real numbers, continuous time evolution 
does not. 

In this paper, we detail and expand the method in 
[l[ which extends the well-known simulation using real 
Hilbert spaces to a more flexible simulation which re- 
quires only local access to ancilla qubits. We also present 
a modification of the Schrodinger equation that allows 
one to view continuous time evolution as a process which 
happens over a real Hilbert space. 

We then present two particular applications where the 
multi-partite simulation answers open questions. Mag- 
niez et al. Q developed an algorithm for testing black 
box quantum devices allowing testing of quantum cir- 
cuits, but the central theorem breaks down for unitaries 
with complex entries. The current work shows that it is 
impossible to test arbitrary unitaries (in their black box 
model) since the real simulation is not unitarily equiva- 
lent to the original system, and yet it produces exactly 
the same measurement results. Hence the theorem can- 
not be strengthened. Also, we answer a question asked by 
Gisin about Bell inequalities. He asked whether any 
Bell inequality could be maximally violated using states 
and measurement operators over real Hilbert spaces. Our 
multi-partite simulation gives an affirmative answer to 
this question. 
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II. SIMULATING QUANTUM SYSTEMS 

A. Simulating complex quantum states and 
operators 

Our goal is to simulate a given quantum system us- 
ing only real numbers in state amplitudes and operator 
matrix entries. This problem can be solved by adding 
one extra qubit and storing the two dimensions of the 
complex plane in the two dimensions of the extra qubit. 
Specifically, we map states as follows (further discussion 
is given in the Appendix): 

J2iax + \x) ^Y,a^\x)\0)+b^\x)\l). (1) 

X X 

In order for interference to happen properly, we need to 
combine the real and imaginary parts correctly. Note 
that the matrix 



XZ = 




applied to the extra qubit corresponds to multiplying the 
original state by i. With this in mind we see that the map 

'^{ax,x'+ibx,x') I a;} (a;' I 1-^ ^ I x) {x \®{ax,x'I+K,x'XZ) 

(2) 

will give operators that map a simulation state of an 
input vector to the simulation state of the correspond- 
ing output vector. Finally, for states | ip) and | (j)) with 
corresponding simulation states | V'') and | (j)') the inner 
products are related by 

(V-V) =Re (^10). (3) 

B. Measurements, Unitaries, and Hamiltonians 

Suppose that we have a positive semi-definite opera- 
tor P with corresponding simulation operator P' given 
as in equation [2l Equation [3] reveals that (V' | -P | V') = 
{ij:' I P' I '0'), hence P' is positive semi-definite. This also 
establishes that POVM outcome probabilities will be the 
same for simulation states and POVM elements as for 
the original states and POVM elements. 

For a unitary U , the corresponding operator U' will 
also be unitary, and any subsequent measurement out- 
comes using a simulated POVM will give also the same 
outcome probabilities as the original states, unitaries and 
POVM. That is to say, we have 

(V- 1 u^PU I V) = (v-' I u'^p'u' 1 1}}') . 



The same substitution will also work for Kraus opera- 
tors of completely positive maps. 

For an analogue of the Schrodinger equation, we ex- 
pand the context a little. Let be a Hermition operator, 
then iH has imaginary eigenvalues and U — exp{iHt) 
will be unitary. Let H' be the simulation operator cor- 
responding to H . Then H' is also Hermitian. Consider 
(/ ® XZ)H': like iH, it will have imaginary eigenvalues 
and hence U' = exp{I XZ)H't will be unitary. By 
checking the Taylor expansion, we can see that U' has 
all real entries, so simulation states will continue to have 
all real amplitudes as they evolve continuously over time 
under this modified Schrodinger equation. 

Combined with the map to simulation states, the sim- 
ulation operators properly evolve the simulation states to 
correspond with the evolution of the original system. In 
particular, the evolved simulation states, combined with 
the simulation measurements, give the same statistics as 
the evolved original states and original measurements. 

Interestingly, in this encoding more operations are pos- 
sible than in the original complex description. In partic- 
ular, complex conjugation can be implemented by ap- 
plying a Z operator to the extra qubit. This allows the 
simulation of arbitrary anti-unitaries as well as partial 
anti-unitaries which act only on a subspace, plus compo- 
sitions of these operators with unitaries. However partial 
transposition, i.e. anti-unitaries acting on only one factor 
of a tensor space, can't be simulated in this way. 

C. Local operations and multipartite systems 

An important aspect of the previous discussion is that 
the simulation of any operation involving an imaginary 
phase involves interacting with the one additional qubit. 
This means that local operations will not in general re- 
main local in the simulation. We can overcome this 
difficulty by storing the additional qubit in a carefully 
constructed two-dimensional subspace stored in several 
qubits. We then add one new qubit to each subsystem in 
the original system and use it to implement local opera- 
tors in the simulation. 

Define the following states on k qubits with y G {0, 1}'^: 

h{y) even 

and 

h{y) odd 

(see the appendix for further discussion) where h{y) is 
the number of ones in y. Interestingly, applying XZ on 
any one of the qubits takes | O) to | l) and | l) to — | O). 
Thus a locally applied XZ on a single qubit has the ef- 
fect of applying a logical XZ on the non-local subspace 
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spanned by | O) and | l). Now instead of using the logical 
XZ in constructing operators, we can choose any conve- 
nient XZ applied to a single qubit out of the fc that 
we have added. Thus local unitaries, Hamiltonians and 
other operators, can be converted into local simulation 
operators. 

III. APPLICATIONS 
A. Self-Testing 

Magniez et. al. (see [2]) prove the following theorem. 

Theorem 1 (Self Testing[§]). Let T he a unitary opera- 
tor with real entries acting on a qubit. Suppose that 

• \ ijj) d A (E) B is a physical state that simulates the 
maximally entangled pair of qubits \ 

• ® /s I V") simulates T ® / | $+) 

• Ga ® Gb I i^) simulates \ $+) 
then Ga is equivalent to T . 

By "equivalence" we mean that there is a unitary that 
maps the physical system to the logical system, plus some 
ancilla, and by "simulates" we means that the measure- 
ment outcome probabilities are identical with respect to 
some specified logical and given physical operators. 

One can imagine that the Hilbert spaces A and B con- 
tain more information than just the logical qubit. Is it 
possible for an adversary to use this additional informa- 
tion to build physical devices that simulate their logical 
description but are not equivalent to it? If T, the logical 
description, has real entries then the theorem tells us that 
this is impossible, but we would like to know whether the 
theorem can be strengthened to include all unitaries. 

Using the simulation developed above an adversary can 
construct the states | O) and | l) on two qubits and store 
them in extra dimensions of A and B. Whether applying 
Ga on A or Gb on B the adversary can perform the real 
simulation of T using local operators that have the form 
required in the theorem. Using this construction all the 
simulation states give the same measurement statistics as 
their ideal counterparts. Note that our real simulation is 
not equivalent to the logical system since inner products 
are not preserved, as seen in equation [3] Thus an adver- 
sary can simulate a complex T in a non-equivalent way 
and there is no stronger theorem (with the same assump- 
tions) for complex T. This is not a counter-example to 
the original theorem because, when the states have real 
amplitudes, the simulation is equal to the logical system 
in a product state with the ancilla systems. 

B. Bell Inequalities 

Gisin asked the question whether Bell inequalities 
could always be maximally violated by states and mea- 



surement operators on a real Hilbert space. Let us con- 
sider the case of two particles and limit ourselves to finite 
dimensional Hilbert spaces. In such a case the Schmidt 
decomposition guarantees that one can choose bases such 
that any state can be written using only real ampli- 
tudes: I V'} = X]j I ''Pj)- Consider now the most 
well known Bell inequality, due to CHSH. Interestingly, 
for any such state, the optimal measurement to that vi- 
olate the CHSH inequality are such that all eigenvectors 
of all observables appearing in the inequality can all be 
written in the Schmidt bases {4>j} and {'0^} using only 
real numbers Q. 

The above observations led some physicists to search 
for inequalities that either require complex numbers (in 
the state and/or eigenvectors of observables) to be vi- 
olated, or for which the use of complex number would 
increase the maximal achievable violation (this would 
be nice: one could experimentally decide whether some 
observed correlations require complex Hilbert spaces or 
whether real ones suffice) . The CHSH inequality requires 
only two measurement on each side. Hence the natural 
next step was to investigate inequalities with 3 measure- 
ments per site, but there is only one such inequality (up 
to symmetries) and again one can chose bases such that 
all amplitudes (state and eigenvectors) are real. Note 
that if one restricts to two dimensional Hilbert spaces 
this has a nice geometrical illustration: measurements 
with real-amplitude eigenvectors are represented on the 
Poincare sphere as vectors all lying on a great circle. Re- 
cently, at last, a Bell inequality for which the optimal set- 
tings do not lie on a circle has been found [11 . Moreover, 
it is elegant: on one side the 3 settings form a orthog- 
onal triedre (like the xyz coordinates) and on the other 
side there are 4 settings on the edges of a tetrahedron. 
Hence, a particular case of the general question outline 
in the previous paragraph is: can this 3x4 setting Bell 
inequality be violated using only real amplitudes by the 
same amount as using complex numbers? 

Navascues et al.Q, using a simulation technique simi- 
lar to the one presented in [ll], and Pal and Vertesi in 
proved that indeed real numbers are sufhcient for 
maximal violation of all bipartite Bell inequalities. The 
present work answers the question in the affirmative for 
inequalities involving any number of parties. However, 
the required dimension increases by a factor of 2 with 
each additional system added. 



IV. CONCLUSION 

We have shown that complex numbers are not required 
in order to describe quantum mechanical systems and 
their evolution, including continuous-time evolution and 
the evolution of multipartite systems. Specifically, there 
is a simulation of such systems using only states and op- 
erators with real coefficients and most importantly this 
can be done with the same assumptions about which sub- 
systems are allowed to interact (that is, with the same 
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local structure). 

We illustrated the non-trivial implications for self- 
testing quantum apparatus and for testing Bell inequali- 
ties. 
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APPENDIX A: WHERE DID THESE STATES 
COME FROM? 

1. Simulation states 



I'/'i) =^k)<»(a. |0}+6, |1)) 

X 



102) = ^\x) ® {-h^\Q) + a^\l)) . 

X 

The first of these is the state given in equation [1] For 
a Hermitian operator P a brief calculation shows that 
Tr(P'p') = {(j)\P' where | (j)) is any normalized lin- 
ear combination of \4>i) and \4>2)- Thus we can replace 
p' with the pure state | ) and get the same outcome 
statistics for any POVM. 

The two dimensional space of pure states that arises by 
applying the simulation map on a pure state corresponds 
to the possible global phases. In fact, by multiplying 
the original state by a phase and applying equation [1] 
one arrives at a superposition of \4>i) and \4>2)- Thus 
arbitrarily choosing | (f)i) instead of another state in the 
subspace is the same as arbitrarily fixing a global phase. 



There is some ambiguity in [T] due to the fact that we 
can multiply a state by a global phase without changing 
the behaviour of the system, but it does result in a differ- 
ent (and possibly orthogonal) simulation state. We can 
resolve this difficulty by looking at the density operator. 
A n X n density operator will be mapped to a 2n x 2n op- 
erator with double the rank and trace 2. To get a density 
operator we thus divide by 2. For a pure state density 
operator, the simulation density operator is rank 2 with 
eigenvalue 1/2 occurring with multiplicity 2. Consider 
the following. For a pure state | tp) = {ax + ihx) \ x) 
(with density operator p = | -0) (-0 |) the simulation den- 
sity operator will be 

p' = i ^ I a;) (x' \®{{axax' + b^hx') I + {bxax' - Uxbx') XZ) 

x.x' 

This can be written as an equal ensemble of the states 



2. Multipartite simulation states 

One way of deriving | O) and | l) is to use the stabi- 
lizer formalism. Note that we want our two dimensional 
subspace to have the property 

-{XZ),(XZ)i\i,)^\i;) 

for all I -0) in the subspace. This just says that applying 
a phase of i independently to two different subsystems 
should be the same as applying a phase of —1. There are 
k — 1 independent operators of this form (for example, 
they can all be written as products of —{XZ)i{XZ)j for 
various j) so there is a two dimensional subspace that 
is stabilized by these operators. The states | O) and | l) 
form a convenient basis for this space. Moreover, this is 
the unique subspace that has the desired property. 



